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We present a systematic study of the effects due to initial condition fluctuations in systems 
formed by heavy-ion collisions using the hydrodynamical simulation code NeXSPheRIO. The study 
was based on a sample of events generated simulating Au+Au collisions at center of mass energy of 
200 GeV per nucleon pair with impact parameter ranging from most central to peripheral collisions. 
The capability of the NeXSPheRIO code to control and save the initial condition (IC) as well 
as the final state particles after the 3D hydrodynamical evolution allows for the investigation of 
the sensitivity of the experimental observables to the characteristics of the early IC. Comparisons 
of results from simulated events generated using fluctuating initial conditions and smooth initial 
condition are presented for the experimental observable elliptic flow parameter (1)2) as a function 
of the transverse momentum, pt, and centrality. We compare 1)2 values estimated using different 
methods, and how each method responds to effects of fluctuations in the initial condition. Finally, 
we quantify the flow fluctuations and compare to the fluctuations of the initial eccentricity of the 
energy density distribution in the transverse plane. 
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I. INTRODUCTION 

In many dynamical systems, the evolution of the sys- 
tem affects its physical characteristics and the informa- 
tion of the early stages is lost, when observing only the 
probes from the final stages. However, by studying how 
much of the initial condition information is retained and 
survives the dynamical evolution it is possible to obtain 
valuable insight not only on the characteristics of the ini- 
tial conditions but also on the interactions that occur in 
the evolution. These conditions can be applied in many 
physical systems such as in the study of the evolution of 
our universe and the anisotropy observed in the CMB if , 
or in the study of the nuclear reactions in high energy 2| . 

In relativistic heavy-ion collisions an important exper- 
imental observable is known as the elliptic flow. Due to 
geometrical anisotropy created in the initial condition by 
the two colliding nuclei and the evolution of the system, a 
space-momentum correlation develops and an anisotropy 
in the azimuthal distribution of the final particles can 
be measured. The elliptic flow or the Fourier expansion 
second harmonic (i^) corresponds to the amplitude of 
this azimuthal anisotropy. Hydrodynamic models pre- 
dict that indeed elliptic flow should be sensitive to the 
eccentricity of the initial conditions Q. Model calcula- 
tions of V2 are also in good agreement with much of the 
experimental data [4-ll2|. Furthermore, two or more par- 
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ticle correlations may carry important information on the 
initial state, hence on the mechanism of QCD dynamics 
at the very early stage of the collision. Recent studies 
[l3l - |2ll ] show that a scenario including fluctuations of the 
initial condition in addition to the geometrical eccentric- 
ity describes better the experimental data than without 
the inclusion of initial state fluctuations. These fluctu- 
ations in the initial condition will affect the eccentricity 
by changing the average value and also by introducing 
an additional fluctuation which will then affect the fi- 
nal observed V2 values. To explore the effects of these 
fluctuations, we need a more precise study, quantifying 
the effects of these initial condition fluctuations to the 
final experimental observables. For this purp ose we have 
used the simulation code NeXSPheRIO [II [HJ that al- 
lows the generation of events using smooth or fluctuating 
initial conditions (IC) and obtain the final particles af- 
ter the hydrodynamical evolution. The procedures used 
in experimental data analysis for V2 determination are 
applied to the simulated data and the effects of the ini- 
tial fluctuations are investigated through the differences 
in the V2 estimates using different methods. NeXSPhe- 
RIO allows for the control and study of both the initial 
condition before the hydrodynamical evolution and the 
final particles after hadronization and freeze-out. Thus, 
it allows for a detailed study, on an event- by- event anal- 
ysis, of the correlation between initial condition param- 
eters such as eccentricity and fluctuations to the final 
state observables such as V2- This simulation code has 
already been extensively tested and has presented rea- 
sonable agreement with experimental data [H, [H, H3] ■ 
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In the next section we describe how the initial energy 
density is obtained in the NcXSPhcRIO code and how 
the eccentricity is calculated. In section HIT1 we describe 
the methods used to estimate the elliptic flow. In section 
IIVI we present our results and, finally, in section [V] we 
summarize our conclusions. 



eccentricity and the participant plane eccentricity as a 
function of the collision impact parameter. In the x-axis 
of Fig. [TJ we also show the equivalent number of charged 
particle density produced at rj = 0, dN c h/drj, which is 
obtained after the full event simulation procedure is com- 
pleted. The top panel shows the event-by-event reaction 



II. INITIAL ENERGY DENSITY 

The NeXSPhcRIO code allows the use of smooth and 
fluctuating initial conditions, which is particularly suited 
and convenient to study the effects of fluctuations in the 
elliptic flow calculations. The initial conditions are gener- 
ated by a microscopic model called NeXuS [22| which pro- 
duces, on an event-by-event basis, detailed space distri- 
butions of the energy-momentum tensor, baryon-number, 
strangeness and charge densities [22j|. The initial condi- 
tion is then used as input to the hydrodynamical model 
and the system evolution is computed up to a given point 
where a decoupling mechanism is applied and final parti- 
cles are produced [22J . Considering the transverse profile 
of the initial energy density distribution around 77 = 0, 
it is possible to calculate the eccentricity of the initial 
geometry. In the situation where there is no fluctuation 
and the colliding nuclei are considered to have smooth 
distributions, the major axis of the almond shape of the 
overlap area at the moment of the collision is perpendicu- 
lar to the plane defined by the impact parameter and the 
beam axis, generally referred to as reaction plane, while 
the minor axis coincides with the direction of the impact 
parameter. 

The eccentricity calculated with respect to this reac- 
tion plane, defined by erp, is given by [2o| : 



b[fm] 
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£rp = 



(1) 



where and a y are the variances of the distribution 
along x and y directions, respectively. Because of fluctu- 
ations, the actual distribution of the hot material created 
by collisions of participants can have the principal axis 
different from those of the smooth initial condition and 
the minor axis deviates from the reaction plane direction 
defined above (25[. The plane that maximizes the eccen- 
tricity determined from the initial energy distribution is 
generally called participant plane 1261 . The participant 
plane eccentricity ep p is given by [25[ : 



£pp 



4< 



(2) 



where a xy is the covariancc of the transverse energy den- 
sity distribution. Note that the reaction plane is deter- 
mined just by the initial geometry of the collision for the 
whole system uniquely, whereas the participant plane de- 
pends on the dynamics of the collisions in each event. In 
Fig. [TJ we show the distributions of the reaction plane 
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FIG. 1. (color online) (a) Reaction plane eccentricity epp 
and (b) participant plane eccentricity epp as a function of 
the number of charged particle produced at r\ = 0. The black 
line in panel (a) represents the mean of the erp distribution, 
which is equal to the epp and epp of the smooth IC case. 
The dashed line in the panel (b) represents the average of the 
epp for the fluctuating IC case. 

plane eccentricity epp and the bottom panel shows the 
participant plane eccentricity epp. By construction, the 
smooth IC reaction plane eccentricity corresponds to the 
mean eccentricity of the fluctuating IC case, represented 
by the solid black line in the panel (a) of Fig. [TJ Also, 
the participant plane eccentricity matches the reaction 
plane eccentricity for the smooth IC case. For the most 
peripheral collisions the reaction plane eccentricity can 
assume negative values due to fluctuation that will force 
the asymmetry plane to be perpendicular to the reaction 
plane. The dashed line in the panel (b) represents the 
mean of the participant plane eccentricity for the fluctu- 
ating IC case only. 



III. ELLIPTIC FLOW DETERMINATION 

The momentum anisotropy of the final particles az- 
imuthal distributions is generally described by a Fourier 
expansion of the azimuthal angle 4> of each particle with 
respect to the reaction plane angle ^pp [27] . The most 
relevant contribution comes from the second coefficient 
of the expansion, V2, which is due to the almond shape 
of the initial overlap area of the incident nuclei in the 
transverse plane. However, the reaction plane angle can- 
not be directly measured experimentally and, therefore, 
several methods and techniques have been developed to 
estimate the vi coefficient from the experimental data 
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[27H30I ] . For this work, we have computed the elliptic 
flow coefficients in several different ways: (a) using the 
true reaction plane (u2{RP}) known in the simulation; 
(b) estimating a reference plane with the event plane 
method (w 2 {EP}) [27|, H3; (c) calculating the participant 
plane from the initial energy density distribution profile 
around r\ = 0; and (d) using the formalism of the cumu- 
lant method (t>2{2} and t>2{4}) [30j . Below we summarize 
how each method works. 

The reaction plane elliptic flow, «2{RP}, is calculated 
with respect to the true reaction plane as: 



w„{RP} = (cos [ntf-Vnp)]), 



(3) 



where <fi is the azimuthal angle of each particle, ^ rp 
is the true reaction plane angle and the angle brackets 
denotes averages first over all selected particles in each 
event, and then over all events. In NeXSPhcRIO, for 
fluctuating IC, an event is first generated by generating 
nucleons inside the colliding nuclei. Thus, before the col- 
lision, the colliding matter already looses its symmetry 
with respect to the reaction plane, as defined by the im- 
pact parameter vector and the beam axis. In this sense, 
the concept of reaction plane looses its physical mean- 
ing for event-by-event calculations. However, we keep 
the calculations of «2{RP} as a reference for comparisons 
with the flow results obtained with other methods. 

The event plane method is used to estimate the true 
reaction plane using the flow anisotropy itself. It is cal- 
culated by defining the flow vector as follows [28[ : 



tan' 1 (Y n /X n ) J X n = J2 l 

uii cos (n<j)i) 

y n = Y l i W i sin ( n( t ) i) 



(4) 



where Wi is a weight and is the azimuthal angle of 
particle i. The flow coefficients are then determined by 



Vn{EP} = 



(cos[n ((/>-#„)]) 
(cos [n (*„ - *rp)]) 



(5) 



where the denominator in the right-hand side of the equa- 
tion [5] is the resolution of the event plane, and the angle 
brackets denote averages taken over all selected particles 
in each event, and then over all events. Experimentally, 
$j;p is not known and thus the resolution must be in- 
ferred in some other way. The technique applied here and 
extensively used in experimental data analysis consists 
in dividing the event being analyzed into two sub-events 
(event a and event b) of same multiplicity and calculate 
the resolution of the sub-event. It has been shown that 
the resolution of the sub-events is given by p7j : 

(cos [n (K - * RP )}) = (cos [n (** - * RP )] ) 



y/(cos[n(n-n)]), (6) 



where and ^ b n are the event plane angles for the sub- 
events a and 6, respectively. From the resolution of the 
sub-events one is able to estimate the resolution of the 



event plane for the full event (see details in Rcfs. |27I . |28| ). 
The event plane angle can deviate from the true reaction 
plane angle due to the fluctuations in the initial condi- 
tions and also due to the resolution caused by the limited 
statistics of the measured particles. Neglecting the angle 
shift due to the statistical resolution, it is expected that 
f bp would follow <f pp, the participant plane azimuthal 
angle. So, in order to check this we also calculate V2 with 
respect to the participant plane angle: 



,{PP} = (cos (n [</>-* PP ])), 



(7) 



where $pp is calculated from the initial energy density 
transverse distribution profile as (25j : 



f pp = tan 



■ni 



4cr. 



■'■;/ 



(8) 
(9) 



A recent study [3l| showed that the event plane angle 
f pp as estimated from the azimuthal particle distribu- 
tion itself seems to be more correlated to the participant 
plane angle ^ pp than to the reaction plane angle VPpp. 
Following this study we have also computed the distribu- 
tions of the correlations between f pp, ^ ep an d ^pp. 
In Fig. [3 we show the distribution of the difference be- 
tween the different angles, calculated for NeXSPhcRIO 
mid-central events. The difference VP pp — f pp (shown 
in red open squares symbols) provides the magnitude of 
the variation in the participant plane direction caused by 
event-by-event fluctuations in the initial energy density 
profile. The distribution of ^ep — ^rp (shown in blue 
open triangles) represents the magnitude of the event 
plane dispersion caused by the IC fluctuations convoluted 
to the effect of resolution caused by the limited number 
of particles used in the ^ ep determination. The distri- 
bution of ^>ep — ^pp (shown in open green circles) is 
formed by the statistical resolution and the width of the 
true correlation between ^ep and 'I'pp, where by true 
correlation we mean the part that comes from the initial 
state fluctuation. For instance, assuming that the par- 
ticipant plane and the event plane are totally correlated 
(as in an event with smooth IC), the distribution of the 
difference between these two angles would be due to the 
statistical resolution only. The width (er) of the distribu- 
tion of the angle difference as a function of the collision 
ccntrality is summarized in Fig. [3] From this plot we can 
see a decrease of the angular resolution for the very pe- 
ripheral collisions and also for the very central collisions. 
The increase of the width for central collision is caused 
by the decrease of the initial eccentricity as well as the 
increase of the lumpiness in the initial energy density dis- 
tribution, therefore, enhancing the effect of fluctuations. 

An alternative method used in experimental data anal- 
ysis to estimate the elliptic flow is through the cumulant 
formalism. The cumulant method calculates the flow co- 
efficients directly from particle correlations, without the 
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FIG. 2. (color online) Correlation between different reference 
plane angles for mid-central NeXSPheRIO events. Red open 
squares (dashed line) show the difference between ^Ppp and 
$iiP, green open circles (dot-dashed line) show the difference 
between ^ep and 'I'pp, and blue open triangles (double-dot- 
dashed line) show the difference between and $bp. 
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FIG. 3. (color online) Width of the correlation between dif- 
ferent reference plane angles as a function of the number of 
charged particles at rj = 0. Red open squares (dashed line) 
show the width of ^pp — ^pp, green open circles (dot-dashed 
line) show the width of ^ ep — ^ pp, and blue open triangles 
(double-dot-dashed line) show the width of ^ep — ^rp- 



explicit need of a reference plane. The prescription used 
in this work is called Q-cumulants or direct cumulants, 
presented in Ref . [3(| . The procedure is divided into two 
parts where in the first part, the reference flow is calcu- 
lated using all particles inside the selected range. In the 
second part, the differential flow (p t dependent) is then 
calculated for the particles of interest. By following the 



notation used in Ref. |3Cj , the flow coefficients for second 
and fourth order cumulants can be written as: 



Reference flow: 



Differential flow: 




v n {2} = V7<2}>_ 



2 -«2» J 



■(10) 



{(2')) 



V«2» ( 
= -[({4'))-2.<(2'))«2))] U 1 -> 

^/-[«4»-2-«2»2] 3 ' 



where the double brackets denote weighted averages of 
2- and 4-particle correlations, first over the particles and 
then over the events. The weights are the total number 
of combinations from two or four particle correlations, 
respectively, and they are used to minimize the effects 
due to multiplicity fluctuations. The advantage of the 
Q-cumulant method is that it is not necessary to per- 
form nested loops to compute all possible combinations 
in multi-particle correlations. Instead, it uses the flow 
vectors to calculate directly the multi-particle cumulants 
(see details in Ref. H(|). The derivation of the expres- 
sions for higher order cumulants is straightforward. We 
present here only up to the fourth order since the sixth 
and higher orders do not seem to differ much from the 
fourth (see for instance Ref. [Hj]). In the next section 
we present the details of the simulated sample used, the 
centrality classes, and the estimates for V2 obtained using 
each one of the methods described above. 



IV. RESULTS 

In order to compare the elliptic flow estimates obtained 
with different techniques, we used a sample of simulated 
events for Au+Au collisions at the center of mass energy 
of 200 GeV per nucleon pair. For the analysis presented 
here we used only the charged particles produced and 
the weak decays have been turned off to minimize the 
non-flow contribution. In our sample each event gen- 
erated has associated a three-dimensional distribution of 
the energy density in the initial overlap region of the inci- 
dent nuclei. The energy density distribution is computed 
from the energy-momentum tensor given by the NeXuS 
code at an initial stage and can fluctuate on an event-by- 
event basis, depending on the nuclear distribution deter- 
mined by the incident nuclei and by the collision impact 
parameter j22|. Alternatively, it is also possible to in- 
put a smooth distribution generated by averaging over 
many events. For the analysis presented in the follow- 
ing we have divided the simulated sample into five event 
centrality classes for both fluctuating and smooth ini- 
tial conditions, as described in table HI A minimum of 
30 thousand events were used for each event centrality 
class, and for both smooth and fluctuating IC. The to- 
tal simulated data sample adds to more than half million 
events. 
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TABLE I. Centrality classes of the NeXSPheRIO events used 
in the study. 



Centrality b ran 


ge (fm) 


(dN ch /dri) 


0-10% 


0.00 


- 4.78 


576.5 ± 0.4 


10-20% 


4.78 


- 6.77 


393.8 ± 0.3 


20-30% 


6.77 


- 8.29 


271.5 ± 0.2 


30-40% 


8.29 


- 9.57 


182.3 ± 0.1 


40-60% 


9.57 - 


- 11.72 


89.9 ± 0.1 



A. Flow comparisons 

We present results obtained for the estimates as 
a function of the transverse momentum, p t , and cen- 
trality of the collision as given by the mean number 
of charged particles produced at 77 = 0. The calcula- 
tions for all methods were performed using charged par- 
ticles within the pseudorapidity window < 1.0 and 
for 0.15 < pt < 2.0 GeV/c. The event plane determi- 
nation was done using particles from 2.5 < | ^7 1 < 4.0, 
therefore, avoiding auto-correlation on -^{EP} calcula- 
tion, and with a requirement of a minimum of 15 charged 
particles within this pseudorapidity region for each event. 
We have also included published results reported by the 
" and PHOBOS m 



STAR [32j, PHENIX [33J and PHOBOS [34J experiments 
for comparison. In the following plots (Figs. 0] to [5]), the 
V2 calculated by the different methods are presented by 
different symbols. Blue open circles for the reaction plane 
V2, red open squares for the V2 calculated through the 
event plane method, gray open diamonds for the partic- 
ipant plane V2, and the green open triangles and orange 
open crosses for the v% calculated using two and four 
particle cumulant methods, respectively. In Figs. [4] and 
[5] we show the results for V2 estimates as a function of 
transverse momentum for both smooth (left panels) and 
fluctuating (right panels) initial conditions, for centrali- 
ties 10-20% and 40-60% respectively. The blue star sym- 
bols are experimental results from the STAR experiment 
[HH and represent V2 estimates obtained with the event 
plane method, and the yellow triangles are results from 
the PHENIX experiment 33j obtained with the second 



order cumulant method. The pt dependence of the V2 
curve from the NeXSPheRIO data generated using fluc- 
tuating IC shows a better agreement to the experimental 
data than the results from the smooth IC. The compar- 
isons presented in Figs. [4] and [5] also show that in the 
case of the smooth IC, the different methods for the V2 
estimate that are based on the measurement of the final 
state particles provide the same results as the t>2 calcu- 
lated from the IC participant plane method. However, in 
the case of fluctuating IC, there is a discrepancy between 
the experimentally measurable V2 methods and the ac- 
tual elliptic flow from the participant plane calculation. 
Moreover, the difference between the results from differ- 
ent V2 methods increases with the transverse momentum. 
This shows clearly how each method is affected by the 
fluctuations in the initial condition. Assuming that the 
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FIG. 4. (color online) Differential V2 as a function of the 
transverse momentum for 10-20% central events. Open sym- 
bols are U2{RP} (blue circles), «2{EP} (red squares), 112 {PP} 
(grey diamonds), V2{2} (green triangles) and «2{4} (orange 
crosses), for smooth (panel (a)) and fluctuating (panel (b)) 
initial condition. Closed blue stars and yellow triangles are 
results from STAR [H and PHENIX [33] experiments, re- 
spectively. The lines are just to guide the eyes. 



participant plane is the plane which defines the direction 
of the elliptic anisotropy created by the initial state and, 
therefore, defines the direction of the V2 in the final state, 
it is possible to take the values of 1>2{PP} as a reference 
for the real (1)2). Comparing to the results obtained with 
methods used experimentally (namely W2{EP}, «2{2} and 
«2{4}), we see f2{EP} above «2{PP} and below v 2 {2}, 
consistent with 



(v 2 ) <v 2 {EP}< v^f) [25]. 

In Fig. [6] we show the dependence of the integrated V2 
with the collision centrality. In the panel (a) we present 
the results obtained with the event plane method, panel 
(b) shows the results obtained with the second order cu- 
mulant, and panel (c) shows the results obtained with the 
fourth order cumulant. Results from the STAR [32[ and 
PHOBOS [3^ experiments are also presented for compar- 
ison. In all three panels we also included the participant 
plane 1/2 as a reference. The behavior of the flow obtained 
with the NeXSPheRIO model is very similar to the exper- 
imental results but the values are systematically lower. 
In Fig. [7] we have combined all the NeXSPheRIO results 
in the same panel and a clear ordering of the values from 
the different methods can be observed. This ordering 
presented by t>2{4}, «2{EP} and voj2} is also observed 
in experimental results (see Ref. (32|). As expect, the 
elliptic flow calculated with respect to the reaction plane 
are lower than all other results, confirming that in a sce- 
nario with lumpy initial energy density distribution, the 
plane defined by the impact parameter vector and the 
beam axis no longer drives the preferred direction of the 
flow. To see how the elliptic flow scales with the eccen- 
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FIG. 5. (color online) Differential «2 as a function of the trans- 
verse momentum for 40-60% peripheral events. Open sym- 
bols are «2{RP} (blue circles), «2{EP} (red squares), ^{PP} 
(grey diamonds), V2{2\ (green triangles) and V2\4} (orange 
crosses), for smooth (panel (a)) and fluctuating (panel (b)) 
initial condition. Closed blue stars are results from STAR 
[3^ | experiment. The lines are just to guide the eyes. 



tricity of the initial state, we plotted the mean integrated 
V2 obtained with the cumulant method over the respec- 
tive eccentricity cumulant moment, calculated using the 
participant plane eccentricity [H, [25|, HH : 



e{2] 
£ {4} 



^ 2 ) 
'PPi 



2(4p) 2 



-PPi 



(12) 
(13) 



as a function of the multiplicity density. Compar- 
ing these ratios with i> 2 {RP} / (s rp) , f2{PP}/(epp) and 
v 2 {EP}/e{2} as shown in Fig. obtained for \rj\ < 1.0 
(fluctuating IC), we find a good agreement among the 
NeXSPheRlO results. The points obtained with differ- 
ent methods seem to fall almost on top of each other, with 
the only exception of ti 2 {2}/e{2}, which is systematically 
higher. Such behavior is being investigated and we have 
already observed that the difference between i> 2 {2}/e{2} 
and the other ratios seems to vanish when increasing 
the pseudorapidity window considered in the calculations 
to |?7| < 6.0. The observable « 2 {4}/e{4} provides quite 
consistent results with ?; 2 {PP}/(£pp), the later one be- 
ing a quantity obtained directly from the anisotropy of 
the IC and, therefore, not accessible experimentally. For 
comparison, we also show results from STAR [3fjj experi- 
ment for v 2 {2}/e{2} and i> 2 {4}/e{4}, with the eccentric- 
ity taken from fKLN-CGC model, for Au+Au collisions 
at the center of mass energy of 200 GeV per nucleon pair. 
The experimental results are always higher than the NeX- 
SPheRlO results, but they present similar separation be- 
tween u 2 {2}/e{2} and u 2 {4}/e{4}. Song et al. [n| have 
recently reported that the inclusion of viscous effects can 
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FIG. 6. (color online) Integrated 1)2 as a function of {dN c h/drj) 
at f] = 0. Panel (a) shows u 2 {EP} for NeXSPheRlO events 
(red open squares) and results from STAR [32| (red solid 
squares) and PHOBOS [34| (green solid circles) experiments; 
panel (b) shows U2\2} for NeXSPheRlO events (green open 
triangles) and from STAR experiment (green solid triangles) 
[321 ] ; and panel (c) shows t>2{4} for NeXSPheRlO events (or- 
ange open crosses) and from STAR experiment (orange solid 
crosses) [32| . The V2 results obtained with respect to the par- 
ticipant plane (grey open diamonds) were also included for 
reference. The lines are just to guide the eyes. 



reduce the baseline of ideal fluid w 2 /e. Moreover, they 
pointed out that a proper event-by-event treatment can 
also affect this baseline. Thus, a study of the observables 
in an event-by-event scenario, as presented here, is also 
important for the comparison between data and viscous 
hydrodynamical models. 



B. Flow fluctuations 

The differences observed for the flow estimates indi- 
cate that each method responds differently to the fluctu- 
ations in the initial condition, which can be used to study 
these fluctuations with the final observables. In particu- 
lar, from the definitions of flow estimates obtained with 
the cumulant formalism, it has been suggested by the au- 
thors of Ref. [1^ that, if only the leading order of a\ 
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FIG. 7. (color online) Mean integrated elliptic flow, (112), 
as a function of {dN c h/dr}) at 77 = 0. Blue open circles are 
«2{RP}, red open squares are «2{EP}, grey open diamonds 
are V2{PP}, green open triangles are «2{2} and orange open 
crosses are t>2{4}. The lines are just to guide the eyes. 
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FIG. 8. (color online) Eccentricity scaled V2 as a function of 
{dN c h/dr)} at r\ — 0. Open symbols are V2{PP}/ (epp) (grey 
diamonds), V2{RP}/(£rp) (blue circles), V2{EP}/e{2} (red 
squares), V2{2} /e{2} (green triangles) and U2{4}/e{4} (or- 
ange crosses), from NeXSPheRIO fluctuating IC. Closed sym- 
bols are results from STAR experiment [35| for H2{2}/e{2} 
(red stars) and f2{4}/£{4} (blue triangles). The error bars 
are the quoted statistical and systematic uncertainties added 
in quadrature. The lines are just to guide the eyes. 



and S arc considered, then: 



V2{2} 2 

V2H} 2 



(v 2 } 2 + 5 + a 2 m 



{V2 



(14) 
(15) 



where S is the non-flow contribution (in general, correla- 
tions other than those related to the reaction plane) , and 
a v . 2 is the elliptic flow fluctuation. The approximation in 
equation 1151 is valid for cr„ 2 <C (^2) and negligible higher 
order moments. We note that this approximation breaks 
down for peripheral and central collisions where the skew- 
ness and kurtosis of the v 2 distribution and terms related 
to {v2)cr V2 will contribute significantly to i>2{4}. From 
equations [14] and [T5] we can extract that the difference 
between «2{2} 2 and «2{4} 2 will provide a quantitative 
measure of the non-flow contribution added to the flow 
fluctuation as given by: 



v 2 {2} 2 - f 2 {4} 2 w 5 + 2a\ 



(16) 



We present this quantity as a function of the transverse 
momentum in Fig. [9l In order to reduce the statisti- 
cal uncertainties in low multiplicity events, the particle 
pscudorapidity acceptance was increased to |t;| < 6.0. 
Results from the smooth IC case, shown in Fig. [9] as the 



20 



15 



x10"' 



CM 
> 
I 

CM 

CM 



10 



5- 



• 0-10% Smooth IC 

O 0-10% Fluctuating IC 

A 40-60% Smooth IC 

A 40-60% Fluctuating IC 



A - 



A 



h 

-2- 



o 



o 
4 



O 

1 

t 



0.0 



0.5 



1.0 
p [GeV/c] 



1.5 



2.0 



FIG. 9. (color online) V2 fluctuations as characterized by 
Eq. [TS]as a function of transverse momentum. Results from 
smooth IC are shown by closed symbols for 0-10% central 
(yellow circles) and 40-60% peripheral (blue triangles) events. 
Results from fluctuating IC are shown by open symbols for 0- 
10% central (circles) and 40-60% peripheral (triangles) events. 

solid yellow circles and the solid blue triangles, for event 
centrality 0-10% and 40-60% respectively, are consistent 
with negligible fluctuation and non-flow contributions in 
the model. Results from the fluctuating IC events are 
presented in Fig. [5] by open circles for the 0-10% cen- 
tral events and open triangles for the 40-60% peripheral 
events. It is clear that in these cases, the quantity defined 
by equation [16] is non-zero, indicating that indeed such 
parameter is sensitive to the fluctuations of the IC. The 
absolute magnitude of the fluctuation is higher in the pe- 
ripheral 40-60% events compared to the central events. 
Fig. ITUI shows the results of ^2{2} 2 — i>2{4} 2 calculated 
for the different event centrality classes. Solid squares 
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FIG. 10. (color online) V2 fluctuations as characterized by Eq. 
1161 as a function of (dN c h/dr)) at 77 = 0. The open diamonds 
are the results for the fluctuating IC case with 77 < 1.0, 
and the squares are the results for the fluctuating IC (open 
squares) and smooth IC (solid squares), both for 77 < 6.0. 
Red solid stars are results from STAR experiment [3q. The 
error bars are the quoted statistical and systematic uncertain- 
ties added in quadrature. 



represent the smooth IC case and open squares the fluc- 
tuating IC case, both for \r]\ < 6.0. The smooth IC case is 
consistent with zero, while the fluctuating IC case shows 
a steady decrease towards the most central collisions. In 
addition, we have also performed the analysis for fluctu- 
ating IC obtained within a tighter pseudorapidity window 
(\r/\ < 1.0), represented by open diamonds in Fig. [10] We 
observe a similar behavior compared to the fluctuating IC 
case for wider r\ window, but with points systematically 
higher. The increase of the values observed for this case 
suggests a dependence of the fluctuations with pseudora- 
pidity window. Published data from STAR experiment 
|35| are also shown for \r\\ < 1.0. These points present 
values always higher than the results from NeXSPheRIO, 
but the behavior is similar. Moreover, the contribution of 
non-flow effects is expected to be greater in experimental 
data. 

Assuming a negligible non-flow component in equation 
[TBI it is possible to take the difference V2{2} 2 — ^2{4} 2 as 
an estimate of the absolute flow fluctuation. An estimate 
for the relative flow fluctuation can also be defined by 
using the equations [141 and [Pol as: 



R v 



1 v 2 {2} 2 - V2H} 2 
« 2 {2} 2 + « 2 {4} 2 ' 



(17) 



In the case that the vi distribution is Gaussian-like, with 
the mean much larger than the width, the quantity R v 
is a resonable approximation for relative flow fluctuation 
<T v 2 / i v 2) [lEHZl- I n Fig. [TJJthe Rv parameter as a func- 
tion of pt is presented for the most central 0-10% and the 



peripheral 40-60% event centrality classes. Opposite to 
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FIG. 11. R v parameter as a function of the transverse momen- 
tum. Open circles represent the 0-10% most central events 
and open triangles are for 40-60% peripheral events. 

what was observed for the absolute flow fluctuations in 
Fig. [9] in this case the central events show higher val- 
ues than the peripheral events. The higher value of R v 
observed in central events compared to the peripheral 
are due to both the increase of the fluctuations in cen- 
tral collisions and also the increase of asymmetry in the 
V2 distributions. The dependence of R v with centrality 
will be further discussed next. In addition, we observe 
negligible dependence with transverse momentum up to 
2 GeV/c for both central and peripheral events. Simi- 
larly to equation [17] it is possible to define an equivalent 
quantity for the eccentricity, R E : 



R, 



/e{2} 2 -e{4} 2 
e{2} 2 + s{4} 2 



(18) 



where £{2} and e{4} are given by the equations [12] and 
I13[ respectively. Since we have access in our simulated 
events to both the IC and also the final state particles, 
we can verify if the quantity R e is a good approximation 
for the relative eccentricity fluctuation, a e /(e). In Fig. 
[12] wc show the R e parameter as defined by equation [T8l 
(blue solid line) and the relative eccentricity fluctuation, 
<J £ /(e) (red dashed line), extracted from the participant 
plane eccentricity distribution of Fig. [TJ panel (b) . We 
have also included for comparison, results for Glauber 
Model (yellow dot-dashed line) and Color Glass Conden- 
sate (CGC) (green double-dot-dashed line) Monte Carlo 
calculations [38j |. In this case, the references reported 
these curves as a function of the number of participants 
in the collision (N var t)- We used the values presented in 
table II of Ref. [39| to convert N part to the equivalent 
(dN c h j drf) at rj = 0, used in this work to summarize the 
collision centrality. This result shows that the observ- 
able R v is a good measure of the relative fluctuations in 
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mid-central collisions, but overestimates the fluctuations 
in central collisions. The higher value of R v in central 
events is mainly due to the contributions of higher mo- 
ments of the £pp distribution such as skewness and kur- 
tosis. The Gaussinan-like condition is no longer satisfied 
and the approximations done in equation [TS] that allows 
for the interpretation that R v is a good measure of the 
flow fluctuation is no longer valid. 





1.0r 




0.9 L 




0.8 L 




0.7 L 








0.6 L 






CC 


0.5 7 




0.4 r 




0.3 r 




0.2 r 








n 1 1 1 1 — i 1 — r 

_ R,. NeXSPheRIO IC 

t pp 

-o e /<£„„> NeXSPheRIO IC 

*PP PP 

. o, /<e ) Glauber MC 

e pati part 

„.a t /<e )CGCMC 

E pan part 




100 200 300 400 500 600 

<dN h /dn> 



FIG. 12. (color online) Comparison between the R E parame- 
ter (blue solid line) and the eccentricity relative fluctuation, 
o>/(e) (red dashed line) calculated with respect to the partic- 
ipant plane in the NeXSPheRIO initial condition, as a func- 
tion of (dN c h/dri) at 77 = 0. It is also shown the relative 
eccentricity fluctuation calculated with Glauber Model (yel- 
low dot-dashed line) and CGC (green double-dot-dashed line) 
Monte Carlo EH- 



As an extension of the work already presented by Hama 
and collaborators [23[ and Sorensen [40]| . we have also 
plotted the R v parameter as given by the equation [T7] as 
a function of the mean number of charged particles pro- 
duced at r} = 0, hence, calculated for the different event 
ccntrality classes. The results are presented in Fig. Q~3] 
Open diamonds were obtained for |rj| < 1.0 and the open 
squares for \r)\ < 6.0. In addition to the R v parameter 
we have also included the R e parameter (blue solid line) 
calculated with respect to the participant plane obtained 
from the IC. Also, results obtained for the dynamic rel- 
ative flow fluctuation, oty n /{v2) (black closed circles), 
reported by the PHOBOS experiment [IH, and the R v 
parameter (red closed stars) reported by the STAR ex- 
periment |35| . are shown for comparison. Therefore, since 
the final state observable that we have is R v , it is impor- 
tant to compare it with the initial state quantity R ei and 
not a £ /{e). The NeXSPheRIO values of for \rj\ < 1.0 
seem to agree very well with the experimental results. 
Even though a perfect agreement between the R e curve 
and the points calculated with the R v parameter is not 
observed, it is remarkable that those quantities, obtained 
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FIG. 13. (color online) Comparison between the R v and R £ 
parameters as a function of (dNch/d-q) at r\ = 0. The symbols 
represent R v for the fluctuating IC case with \r/\ < 6.0 (open 
squares) and \t]\ < 1.0 (open diamonds). The blue solid line is 
the R e parameter calculated with respect to the participant 
plane in the NeXSPheRIO IC. Closed symbols are published 
results from PHOBOS [H (black circles) and STAR [U (red 
stars) experiments. The shaded band represent the errors 
quoted by PHOBOS and the error bars in STAR data are the 
statistical and systematic uncertainties added in quadrature. 



using the properties of the very beginning and the very 
final stage of the system evolution, still yield such similar 
values. Thus, we can conclude that the experimentally 
observable parameter R v is a good estimate of the fluctu- 
ations in the initial state of the collisions, for mid-central 
events, which cannot be probed directly. Moreover, we 
can conclude that the R v parameter does not get affected 
by the hydrodynamic evolution and the freeze-out pro- 
cess. 



V. SUMMARY 

In this work we have analyzed a large amount of sim- 
ulated events produced with the NeXSPheRIO code and 
both the IC and the final state particles were saved and 
analyzed on an event- by- event basis. Smooth and fluctu- 
ating IC were used to generate the events for five differ- 
ent centrality classes. This allowed us to study the effects 
of fluctuations in the initial energy density distribution 
through the 1)2 estimates from the final particle azimuthal 
distributions. We showed for the smooth IC cases that 
the methods generally used in experimental data anal- 
ysis to estimate i>2, namely ^{EP}, f2{2} and «2{4}, 
produce consistent results with flow estimates that uses 
informations from the initial condition, namely ^{RP} 
and 1>2{PP}- Also, for the fluctuating IC we found that 
the methods start to deviate from each other when going 
to higher transverse momentum, giving rise to a system- 
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atic ordering as a function of centrality. The discrepancy 
between the different V2 calculation methods were ex- 
ploited to be used as a measurement of the IC degree 
of fluctuation. Our results show that the magnitude of 
the non-flow plus the flow fluctuations increases with p t 
for both central and peripheral events, being more pro- 
nounced for the later one. On the other hand, the rela- 
tive fluctuations, as defined by the R v parameter, shows 
no dependence with transverse momentum. We also ob- 
served a dependence of the quantity defined in equation 
1161 with the collision centrality and the pseudo-rapidity 
window of the particles used in the analysis. Although 
the behavior observed for R v and R e as a function of the 



mean number of charged particles at 77 = do not com- 
pletely agree, the values obtained are remarkably similar, 
which indicates that indeed the measurements of the final 
state flow fluctuations may provide important informa- 
tion on the initial state fluctuations. 
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